BCS superconductivity of Dirac electrons in graphene layers 
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Possible superconductivity of electrons with the Dirac spectrum is analyzed using the BCS model. 
We calculate the critical temperature, the superconducting energy gap, and supercurrent as functions 
of the doping level and of the pairing interaction strength. Zero doping is characterized by existence 
of the quantum critical point such that the critical temperature vanishes below some finite value of 
the interaction strength. However, the critical temperature remains finite for any nonzero electron 
or hole doping level when the Fermi energy is shifted away from the Dirac point of the normal-state 
electron spectrum. We analyze the behavior of the characteristic length scales, i.e., the London 
penetration depth and the coherence length, which determine the critical magnetic fields. 
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Graphite attracts attention of experimentalists and 
theorists for a long time. The interest is explained by 
unusual properties of this quasi-two dimensional mate- 
rial, which are mostly related with the existence of the 
Dirac or conic point in the electronic spectrum (see Fig. 
[T]). Though the theory has predicted the existence of 
such point in graphite many decades ago[l], only re- 
cently the experimental evidences of its existence were 
received: first in graphite [2], which is believed to be a 
stack of weakly coupled atomic layers, and soon after it 
in graphene [1, 0|. The latter discovery has triggered an 
avalanche of experimental and theoretical works. More- 
over, graphene can display unusual properties as a part 
of normal-superconducting hybrid structures: For exam- 
ple, the Andreev reflection has been predicted to have 
new features not characteristic for typical contacts [5]. 

Thorough investigation of graphite has revealed also 
evidences of intrinsic superconductivity in doped sam- 
ples (see Refs. [6|, l7j and references therein). Various 
mechanisms of superconductivity in graphene have been 
considered theoretically. Phonon and plasmon mediated 
mechanisms were discussed in Ref. [8] whereas a res- 
onating valence bond model was proposed in Ref. [9[. 
The Cooper pairing in the undoped graphene may expe- 
rience problems because the Fermi surface shrinks near 
the Dirac point and reduces to zero the number of states 
at the Fermi energy. Indeed, it was shown within the BCS 
model [To[ that the superconducting transition in the un- 
doped graphene possesses a quantum critical point at a 
finite interaction strength below which the critical tem- 
perature vanishes. However, one would expect that the 
electrons in graphene may become unstable towards for- 
mation of Cooper pairs for any finite pairing interaction 
if doping shifts the Fermi level away from the Dirac point 
because the behavior of electrons in the latter case bears 
more resemblance to that in usual metals. This idea has 
been qualitatively discussed in Refs. [8|, LlJj and verified 
within the resonating valence bond model in Ref. 



The aforementioned investigations of superconductiv- 
ity in graphene or graphite (except for Ref. [Io|, where 
only the undoped case was considered) were done taking 
into account the specific details of each particular pairing 
mechanism. However, it would be worthwhile to perform 
the analysis in a more general form independent of a par- 
ticular nature of the pairing mechanism. In the present 
Letter we apply the standard s- wave BCS model for the 
Dirac spectrum of electrons with a minimum number of 
parameters characterizing the pairing interaction, i.e., its 
intensity and the range of interaction in the momentum 
space. The values of these two parameters may vary de- 
pending on the mechanism. Such approach inevitably 
ignores some details and thus is less accurate. However, 
we hope that the loss of accuracy is compensated by a 
more general and transparent picture of the most essen- 
tial features of the Cooper pairing in the systems with 
the Dirac spectrum. 

In what follows we calculate the critical temperature, 
the superconducting energy gap, and the supercurrent as 
functions of the doping level and the pairing interaction 
strength. Without doping the critical temperature van- 
ishes below some finite value of the interaction strength. 
However, the critical temperature is nonzero for any 
nonzero electron or hole doping level when the Fermi en- 
ergy is shifted from the Dirac point of the normal-state 
electronic spectrum. This provides the quantitative ba- 
sis for the earlier conjectures of Refs. M, [ll[ and agrees 
qualitatively with the results of Ref. [9] for resonating 
valence bond model. Moreover, by analyzing the effect 
of the Dirac point on the supercurrent we demonstrate a 
novel feature that, as distinct from the usual supercon- 
ductors, the supercurrent density is not proportional to 
the total number of electrons but is drastically decreased 
due to the presence of the Dirac point. Finally, we esti- 
mate characteristic length scales (penetration depth and 
coherence length), relevant for determination of the crit- 
ical magnetic fields. 
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FIG. 1: (Color online) Conical energy spectrum, (a) Un- 
doped, (b) electron-doped, and (c) hole-doped spectrum 



Consideration of a two-dimensional model requires a 
few comments concerning the applicability of the mean- 
field approach. It is well known that the superfluid tran- 
sition in a two-dimensional system occurs in the form of 
the Berezinskii-Kosterlitz-Thouless transition at a tem- 
perature lower than the mean-field transition tempera- 
ture. Therefore our calculations provide the upper bound 
on the critical temperature and give a good estimate for 
the temperature scale of the transition Moreover, 
the applicability of the mean-field approach improves for 
graphite where a nonzero interplanar coupling, however 
small, does always exist. 

Spectrum. - We assume the energy spectrum in 
graphene in the form 
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The upper or lower sign refers to the conduction or va- 
lence band, respectively; Ep is the Fermi energy without 
doping when the Dirac point lies at the Fermi level. If 
the Fermi energy is shifted by some amount \i due to 
doping, Ep = Efo + /i, (see Fig. [I]) the energy measured 
from the Fermi level is 



Ep = ±vp — \i . 



The group velocity is d^ p /dp = ±vn where n 
electron doping, \i > 0, we have for £ p < 



P : 



p/p. For 



(i) 



Similar relations take place in the case of hole doping, 
fj, = — |/i|, as well. 

BCS gap equation. - We use the standard BCS theory 
and assume an s-wave pairing interaction V p = — \V P \, 
where \V P \ ~ \V\a 2 is the Fourier transform of the pairing 
potential, V is the energy amplitude, and a is the range of 
interaction. We do not consider here the nature of pai ring; 
interaction, but refer the reader to Refs. 
where various possible mechanisms are discussed. The 
coupling constant A is introduced through the equation 
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Here £ TO determines the interval where attractive inter- 
action is present. In what follows we consider the case of 
low doping when < £ TO . 

The BCS gap equation in a spatially uniform case is 
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where the energy of excitations is E p = y £ 2 + A 2 , the 

phase volume is d 2 p = pdp d(j), where <p is the azimuthal 
angle of n, and n(E p ) is the equilibrium Fermi distribu- 
tion of quasiparticles with energies E p . For zero temper- 
ature 1 - 2n(E) = sign(£). With help of Eq. {J), the 
BCS gap equation becomes 
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for both electron and hole doping. 
For zero doping [i = we have 

Ao = £m(A 2 - 1)/2A . 



(3) 
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Nonzero A is possible only for the strong-coupling limit 
A > 1 [10|. However, Eq. (|3]) shows that, for a finite 
doping, a finite Ao exists even in the weak coupling limit 
A < 1. In the case of low doping level when Ao is small, 
Eq. (|3]) gives the gap in a BCS form 



A = 2|/x| 
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with the prefactor determined by the doping level \fi\ 
rather than by the range of interaction. 

Temperature dependence. - For a finite temperature 
Eq. ([2]) yields the gap equation 



% = 2Tln 
A 



coshjyg + A 2 /2T) 
cosh( v /M 2 + A 2 /2T) 
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For T^Owe return to Eq. Equation © leads to 
the equation for the critical temperature 



$(W2T c ;A)=F(|/i|/2T c ) 



(7) 



where 



<E> (y; A) = A 1 y — In (cosh?/) 

nx 

F{x) = x I (x') -1 tanha/ dx r — ln(coshx) 
Jo 

where F(x) > 0. The critical temperature found from 
Eq. © is plotted in Fig. [2j 
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FIG. 2: (Color online) Normalized critical temperature t — 
2T c /^m as a function of the interaction constant A for various 
doping levels |/x|/£ m . The quantum critical point is at A = 1 
and — 0. 



For ji = the critical temperature satisfies 
$(£ m /2T c ;A) =0, i.e., 



£ m /A = 2T c ln [cosh(£ m /2T c ; 



(8) 



This equation has a solution only for interaction strength 
above the quantum critical point, A > 1 (see Fig. [2]). If 
A — > 1 we have T c = £ m (A — l)/21n2 which vanishes 
at A = 1. Comparing this with Eq. (j4]) we find Ao = 
T c 21n2. In the other limit A > 1 we find T c = £ m A/4 
and Ao = 2T C . These results agree with Ref.jlQj where 
only undoped case was considered. 

However, for any low but finite doping level the critical 
temperature is finite. Consider weak coupling limit A <C 
1 where we expect T c <C fi. Indeed, the l.h.s. of Eq. ([7j) 
is $(y; A) = 2/(A _1 - 1) + In 2 already for T c < f m . On 
the other hand, for x ^> 1 the r.h.s. of Eq. ([7j) is 
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where 7 = e c = 1.78 and C = 0.5772 is the Euler con- 
stant. This yields 

rr, 2|/i| 7 



■ exp 
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resulting in the BCS relation A = (tt/j)T c = 1.76T C . 

Consider the vicinity of the quantum critical point \i = 
and A = 1. On the weak coupling side A < 1, the critical 
temperature is given by Eq. ([TO]) which is exact provided 
T c <C i.e., for 1 - A > |/i|/£ m . For |/i|/£ m ~ 1 
and A —> 1, Eq. ([TO]) works also reasonably well. For 
example, Eq. (TOD gives T c w 0.42|/i| for A = 1. This 
can be compared to the exact value for A = 1 which is 
found from the condition F(\/j,\/2T c ) = In 2 resulting in 
T c « 0.40|/i|. In the limit |W « T c « ( m which is more 
appropriate on the strong-coupling side of the quantum 
critical point 



MA - 1) + N /e(A-l) 2 + M 2 21n2 



This holds for |/i|/£ m <C A — 1 <C 1, but also matches with 
the exact T c by the order of magnitude when A — > 1. 

Therefore we come to the conclusion that a finite T c 
does always exist for a finite fi. If A > 1, the critical 
temperature is close to that determined by Eq. (|8|) as 
long as ji <C £ m . If A < 1 we essentially have Eq. ([TO]) . 

Supercurrent. - Let us assume a homogeneous flow 
of the condensate: A = |A|e* ksF , where k s = Vx is 
a constant gradient of the order-parameter phase. In 
the presence of magnetic field, k s = Vx — (2e/hc)A. 
Consider the state described by the particle-like and hole- 
like Bogoliubov-de Gennes wave functions 

u(r) = u p e ip +- r/h , v(r) = v p e ip ~- r/h , 

where p± = p ± Kk s /2, 

e p = e d + em , E M = yJ$ + \Kp, 

Ejj = (d£ >p /dp)h'k s /2 is the Doppler energy, and 

Up = + w4 0) ) 1/2 > ^ = - w4 0) ) 1/2 

are the coherence factors. The standard expression for 
the current is 
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Expanding Eq. (fTT|) for small £x> <C A, T and mak- 
ing shift in the momentum variable we find for the two- 
dimensional current density in the linear response regime 
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4ir 2 H 2 dp 
This yields the current 

j = (eA/47rfi)k, 



where A is the characteristic energy. For zero tempera- 
ture we have 



A = 2IAI 



41n2 



In contrast to the usual superconductors the supercurrent 
density is not proportional to the total electron density, 
being drastically affected by the presence of the Dirac 
point. In particular, for weak coupling limit, |A| <C /i, 
the current j = e/ik s /47rft is proportional to \fi\ oc y/n, 
where n is the density of free carriers provided by dop- 
ing. Near the quantum critical point when T <C £ m , the 
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current is determined by the superconducting gap itself. 
Indeed, for zero doping, Eq. ([TT]) yields 

A= lAltanh^ . 
i i 2T 

For low temperatures, T <C |A|, we have j = e\ A|k s /47r/L 
Close to T c , where |A| <C T c , the current assumes the 
Ginzburg-Landau form j = e\ A| 2 k s /87rfiT c . 

Characteristic lengths and critical fields. - As usual 
the critical fields are determined by two spatial scales: 
the coherence length £o and the London penetration 
depth Xl. At zero temperature the London penetration 
length for a graphene layer with thickness d is 

A -2 _ 2e 2 A 

h 2 c 2 d ' 

It diverges near the quantum critical point A — > 1, 
/i — > 0. For the undoped case the London length is 
Xl = ($>o/ 7r )'\/d/2A, where $o = Trhc/e is the magnetic- 
flux quantum. Close to the critical temperature the 
London length, Xl = ($>o/7rA)y/T c d, is inversely pro- 
portional \JT C — T as in conventional superconductors. 
The coherence length has a standard form: £o ~ hvp/ A. 
Thus the Ginzburg-Landau parameter k = Al/^o, which 
characterizes the type of superconductivity, does not de- 
pend on the temperature near T c as is the case in con- 
ventional superconductors: 

Xl c l T c d 
£o v F ]l e 2 

For typical values vp = 10 8 cm/s, d = 10 -7 cm, and 
for T c ~ 1 K, the Ginzburg-Landau parameter is on the 
border between the two types, k ~ 1. Therefore, close to 
the quantum critical point where T c — > the supercon- 
ductivity definitely becomes of type I. 

To summarize, we have calculated the critical temper- 
ature, the superconducting gap, and the super current 
as functions of the doping level and of the interaction 
strength for an s-wave pairing within the BCS model. 
The superconducting transition in the undoped graphene 



has a quantum critical point with respect to the interac- 
tion strength, which disappears for any finite doping level 
such that a finite critical temperature exists for any weak 
pairing interaction. The amplitude of the supercurrent 
is drastically affected by the presence of the Dirac point, 
which leads to non-trivial behavior of the characteristic 
length scales (penetration depth and coherence length) 
determining critical magnetic fields. 

We thank Y. Kopelevich and V. Eltsov for stimulating 
discussions. This work was supported by the Forscheimer 
Foundation of the Hebrew University of Jerusalem, by 
the Academy of Finland (grant 213496, Finnish Pro- 
gramme for Centers of Excellence in Research 2002- 
2007/2006-2011), by the ULTI program under EU con- 
tract RITA-CT-2003-505313, and by the Russian Foun- 
dation for Basic Research grant 06-02-16002. 



[1] P. R. Wallace, Phys. Rev. 71, 622 (1947). 

[2] I. A. Luk'yanchuk and Y. Kopelevich, Phys. Rev. Lett. 

93, 166402 (2004). 
[3] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 

A. A. Firsov, Nature 438, 197 (2005). 
[4] Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 

438, 201 (2005). 
[5] C.W.J. Beenakker, Phys. Rev. Lett. 97, 067007 (2006). 
[6] R. R. da Silva, J. H. S. Torres, and Y. Kopelevich, Phys. 

Rev. Lett. 87, 147001 (2001). 
[7] Y. Kopelevich, S. Moehlecke, and R. R. da Silva, in Car- 
bon Based Magnetism, edited by T. Makarova and F. 

Palacio (Elsevier Science, 2006), Chap. 18. 
[8] B. Uchoa and A. H. Castro Neto, Phys. Rev. Lett. 98, 

146801 (2007). 

[9] A. M. Black-Schaffer and S. Doniach, Phys. Rev. B, 75, 
134512 (2007). 

[10] E. C. Marino, and Lizardo H.C.M. Nunes, Nuclear 
Physics B, 741, 404 (2006); Physica C 460-462, 1101 
(2007); Nuclear Physics B, 769, 275 (2007). 

[11] E. Zhao and A. Paramekanti, Phys. Rev. Lett. 97, 230404 
(2006). 

[12] K. Sasaki, J. Jiang, R. Saito, S. Onari, and Y. 
Tanaka, J. Phy s. Soc. Jpn. 76, 033702 (2007); 
|arXiv:cond-mat /06 1 1452 



